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K + final states are studied using a sample of 14 x 10 6 ij{2S) decays 
collected with the Beijing Spectrometer (BESII) at the Beijing Electron-Positron 
Collider. The branching fractions of ij(2S) decays to A" + A"~7r + 7r“7r°, u)K + K~, 
w/o(1710), AT*(892) 0 AT _ 7r + 7r° + c.c., if*(892)++ c.c., AT* (892)+AT" p° + 
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c.c. and AT* (892)° A" - p + + c.c. are determined. The first two agree with previous 
measurements, and the last five are first measurements. 

PACS numbers: 13.25.Gv, 14.40.Cs, 13.40.Hq 


I. INTRODUCTION 


From perturbative QCD (pQCD), it is expected that both J/0 and 0(25) de¬ 
caying into light hadrons are dominated by the annihilation of cc into three gluons, 
with widths proportional to the square of the wave function at the origin |T(0)| 2 [lj. 
This yields the pQCD “12%” rule: 

Bp(2S)^h -®V(2S)—>e+e“ 


Qh ~ 


Bj/ip—*h >e+e“ 

The violation of the above rule was first observed in the pi r and K* + K~ +c.c. decay 
modes by Mark-II j^j]. Following the scenario proposed in Ref. [iil, that the small 
0(25) —> pn branching fraction is due to the cancelation of the S- and D-wave 
matrix elements in 0(25) decays, it was suggested that all 0(25) decay channels 
should be affected by the same S- and D-wave mixing scheme, and thus all ratios 
of branching fractions of '0(25') and J/0 decays into the same final state could 
have values different from 12%, expected between pure IS and 2S states 4], The 
mixing scenario also predicts 0(3770) decay branching fractions since the 0(3770) 
is a mixture of S- and D-wave charmonia, as well. Many channels of J/0, 0(25) 
decays, and 0(3770) decays should be measured to test this scenario. 

In this paper, we report measurements of K + K~ tt + tt~-k 0 final states, as well as 
some intermediate states that decay to the same final states. The data samples used 
for this analysis consist of 14.0 x 10 6 (1 ± 4%) 0(25) events taken at a/s = 3.686 
GeV [5] and 6.42(1 ± 4%)pb _1 of continuum data at a/s = 3.65 GeV 


12 %. 
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II. BESII DETECTOR 

BESII is a large solid-angle magnetic spectrometer which is described in de¬ 
tail in Ref. 7]. The momentum of charged particles is determined by a 40- 
layer cylindrical main drift chamber (MDC) which has a momentum resolution of 
(Xp/p=1.78%\/l + p 2 (p in GeV/c). Particle identification (PID) is accomplished us¬ 
ing specific ionization ( dE/dx ) measurements in the drift chamber and time-of-flight 
(TOF) information in a barrel-like array of 48 scintillation counters. The dE/dx 
resolution is (JdE/dx — 8.0%; the TOF resolution for Bliabha events is <Jtof — 180 ps. 
Radially outside of the time-of-flight counters is a 12-radiation-length barrel shower 
counter (BSC) comprised of gas tubes interleaved with lead sheets. The BSC mea¬ 
sures the energy and direction of photons with resolutions of cr#/ E ~ 21 %j\/~E (E 
in GeV), 07 = 7.9 rnrad, and a z = 2.3 cm. The iron flux return of the magnet 
is instrumented with three double layers of proportional counters that are used to 
identify muons. 

A GEANT3 based Monte Carlo (MC) simulation package jsj], which simulates 
the detector response, including interactions of secondary particles in the detector 
material, is used to determine detection efficiencies and mass resolutions, as well 
as to optimize selection criteria and estimate backgrounds. Reasonable agreement 
between data and MC simulation is observed in various channels tested, including 
e + e _ —► ( 7 )e + e _ , e + e _ —> ( , y)p + p~, J/%b —*► pp, J/i/i —> pi r, and , 0(2S') —> 7 r + 7 r 
J/iji -»• l + l~. 

III. EVENT SELECTION 

The K + K~ir + n~ir° final states are reconstructed with four charged tracks and 
two photons. 


A. Photon and charged particle identification 

A neutral cluster is considered to be a good photon candidate if the following 
requirements are satisfied: it is located within the BSC fiducial region, the energy 
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deposited in the BSC is greater than 50 MeV, the first hit appears in the first six- 
radiation lengths, the angle between the cluster development direction in the BSC 
and the photon emission direction is less than 37°, and the angle between the cluster 
and the nearest charged particle is greater than 15°. 

Each charged track is required to be well fit by a three dimensional helix, to 
originate from the interaction region, V xy = ^jV x + V y < 2.0 cm and \V Z \ < 20 
cm, and to have a polar angle | cos#| < 0.8. Here V x , V y , and V z are the x, y, z 
coordinates of the point of closest approach of the track to the beam axis. The 
TOF and dE/dx measurements for each charged track are used to calculate Xpid (®) 
values and the corresponding confidence levels Probpip>(i) for the hypotheses that 
a track is a pion, kaon, or proton, where i (i — ir/K/p) is the particle type. 

B. Selection criteria 

For the final states of interest, the candidate events are required to satisfy the 
following selection criteria: 

1. The number of charged tracks in the MDC equals four with net charge zero. 

2. The number of good photon candidates equals two or three. 

3. For each charged track, the particle identification confidence level for a candi¬ 
date particle assignment is required to be greater than 1%. 

4. The angle between two photons satisfies 0 77 > 6° to remove the background 
from split-off fake photons. 

5. To reduce contamination from il>(2S) —> with rj —> 7r + 7r _ 7r° and J/ip —> 

y + , N+ lt + N^ lt < 4 is used, where N^ lt is the number of muon identification 
layers with matched hits for the higher momentum positive charged track and 
ranges from 0 to 3, indicating not a muon (0) or a weakly (1), moderately (2), 
or strongly (3) identified muon j^|, and N^ lt is the corresponding number for 
the higher momentum negative charged track. 



6 


6 . To reject background from 0( 2S ) —> 7 r + 7 r J/0, |Mrecoil ~ 3.097| > 0.05 
GeV/c 2 is used, where is the mass recoiling against the 7 r + 7 r - pair. 

To improve track momentum resolution and reduce background, four constraint 
kinematic fits imposing energy and momentum conservation are performed. We loop 
over all combinations of charged tracks and good photons and select the one with the 
minimum Xcom f° r the assignment 0(25) —> r y r yn + n~ K + K~ , where the combined 
X 2 , Xcom; i s defined as the sum of the x 2 values of the kinematic fit, Xkinei an d those 
from each of the four particle identification assignments: xlom = J2Xpid(^) + Xkine ■ 

i 

We require Xcom( 77 ^ K+K ) < 7T+7T 7T+7T ), K+K ) < 

xloJ'nK+K-K+K-), and the confidence level of the kinematic fit to be greater 
than 1 %. 

After applying the above selection criteria, there is still background remaining 
from 0(25) —> KsK^^n+n 0 , Ks —■► 7 r + 7 T _ . This is removed by requiring rn 7r + 1 T - > 
0.51 GeV/c 2 or m n + n ^ < 0.45 GeV/c 2 , rriK+TT- > 0.51 GeV/c 2 or ttik+tt- < 0.47 
GeV/c 2 , and m K - 7T + > 0.51 GeV/c 2 or rriK-n+ < 0.47 GeV/c 2 , where in calculating 
rriKit, the 7 r mass is used for the K ± track. 


IV. DATA ANALYSIS 

The 77 invariant mass distribution for events that survive the selection criteria is 
shown in Fig. Qa), where a clear 7T° signal can be seen. By fitting this distribution 
with a 7T° signal shape, obtained from MC simulation, and a 3rd order background 
polynomial, 698 ± 41 events are obtained in the a/s = 3.686 GeV data sample. 
From exclusive Monte Carlo simulation, we determine that the main background 
channels are from ^(2 S) —> 7 r° 7 r° J/-/, J/-0 —> K+ K~tt+tt~ , il>(2S) — > 7 XcJiXcJ ~ 
K+K~n+7T~, and - 0 ( 2 S) — > 7 Xcj,Xcj K+K~n+n~n°. However, the 77 invariant 

mass distributions from all these background channels do not have a peak at m 77 = 
m n 0 , and, therefore, they will not contribute to the number of fitted events. Applying 
the same selection criteria to the a/s = 3.65 GeV data, we obtain the 77 invariant 
mass distribution shown in Fig. mb). Fitting Fig. □;b) in a similar way as Fig. da) 
yields 35 ±7 events. The efficiencies are (3.68 ±0.05)% for '0(25) —■> K+ K~tt+tt~tt 0 
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and (2.59±0.05)% for continuum e + e - —> K + K~tt + 7t~tt 0 : where the difference is due 
to the initial state radiation for continuum data, ffere in calculating efficiencies, we 
have considered the contributions from important intermediate states (shown later 
in the paper) in ^(2S) —> K + K~n + n~ / K°. 




-fl invariant mass (GeV/c 2 ) 


FIG. 1: The 77 invariant mass distributions of (a) selected —» K + K~TT + TT~'yy 

candidate events and (b) selected e + e _ —* I\ + K~ir + ir~'y'y candidate events at y% = 3.65 
GeV. The squares with error bars are data, the histograms are the fit, and the dashed 
curves are background shapes from the fit. 

I 11 the following analysis, we study the 7 r + 7 r _ 7 r°, Ii + K~, K + n~ + c.c., K + tt° + 
c.c. and 7 r + 7 r~ invariant mass spectra to look for possible intermediate resonance 
states. Fig. Ha) shows the 7r + 7r 7 t° invariant mass distribution for K + K 7r + 7r 77 
events after requiring |m 77 — 0.1351 < 0.03 GeV/c 2 (to increase the efficiency, the 
requirements on m n + n - and m Kn used for the K s veto are removed); a clear to 
signal is observed. A fit with an to signal shape obtained from MC simulation and a 
polynomial background gives 78 ± 11 signal events with a statistical significance of 
8.0a. The detection efficiency for this decay mode is (2.66 ± 0.04)%, where we have 
considered the contributions from important intermediate states. Fig. m shows 
the n + 7r~7r° invariant mass distribution for y/s = 3.65 GeV data, where we obtain 
Olg 3 toK + K~ events at the 68.3% confidence level with a similar fit to the 7 r + 7 r _ 7 r° 
invariant mass spectrum. 


























FIG. 2: The 7r + 7r _ 7r° invariant mass distributions of (a) selected ip(2S) —► K + 
candidate events and (b) selected e + e _ —> K + I \candidate events at yfs = 3.65 
GeV. The squares with error bars are data, the histograms are the fit, and the dashed 
curves are background shapes from the fit. 

The Dalitz plot of events satisfying \m- K +Tr--K° ~ m ui\ < 0.04 GeV/c 2 is shown 
in Fig. Ha). There are vertical and horizontal directions bands corresponding to 
/li(1270) ± —> and a diagonal band corresponding to a K + K~ resonance. 

After requiring m W K > 1-5 GeV/c 2 to remove the /l 1 (1270) ± , the K + K~ invariant 
mass distribution is shown in Fig. Hb). There is a cluster of events near 1.7 GeV/c 2 , 
and there is no similar structure in u sideband events (0.06 < m 7V + n - 7 r o — m u < 
0.1 GeV/c 2 or — 0.1 < m 77 + 7 r-^o — m u < —0.06 GeV/c 2 ), as shown in the hatched 
histogram. Assuming the peak is fo(1710), a fit with the mass fixed at 1.714 GeV/c 2 
and width fixed at 140 MeV/c 2 |Hj| gives 18.9±6.2 events. The statistical significance 
for ij}(2S) —> u;/ 0 (1710) is 3.7cr, and the detection efficiency for this decay mode is 
(2.62 ± 0.06)%. No signal for e + e“ —> o;/ 0 (1710) is observed at y/s = 3.65 GeV. 

After requiring \m n + 7T - 7r o — 0.783| > 0.04 GeV/c 2 to remove ip(2S) —> uK + K~ 
events, the K ± n T mass distributions of K + K~'k + tt~7t° candidates are shown in 
Fig.QJfor (a) y/s = 3.686 GeV data and (b) y/s = 3.65 GeV data. By fitting the Ktt 
invariant mass spectrum with the Monte Carlo determined shape for signal, the Iin 
mass distribution of 7r° sideband events (0.18 < m 77 < 0.21 GeV/c 2 or 0.06 < m 77 < 
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FIG. 3: (a) The Dalitz plot of selected i/j(2S) —> ujK + K~ candidates, and (b) the K + K~ 
invariant mass distribution for candidate ip(2S) —> u)K + K~ events after the A'i(1270) 
veto (■ m W K > 1-5 GeV/c 2 ). The hatched histogram is from u> sidebands, and the curves in 
(b) show the fit described in the text. 

0.09 GeV/c 2 ) to describe the peaking background (mainly y K* (892)° K~n + + c.c.), 
and a Legendre polynomial for other backgrounds, as shown in Fig. Ha), 281 ± 
30 events at a/s = 3.686 GeV for K*(892) 0 K~tt + 7t° + c.c. are obtained. For the 
continuum data, the peaking background is negligible, and the fit is performed with 
a signal shape and a Legendre background polynomial and yields 15 ± 5 events, as 
shown in Fig. lib)- Similarly, the K ± 7T° mass distributions are shown in Fig. H and 
150±26 events at a/s = 3.686 GeV and 6±5 events at a/s = 3.65 GeV are obtained 
for K*{892) + K~ti + ti~ ± c.c., as shown in Figs. H a ) an d (b). After subtracting the 
continuum contributions, we obtain 238±34 events for ip(2S) —» K* (892)° K~n + 7T°± 
c.c. with a detection efficiency of (3.00 ± 0.06)% and 133 ± 30 events for ip(2S) — ■> 
K* (892) + K~n + 7r~ + c.c. with a detection efficiency of (3.02 ± 0.06)%. 

The n + 7i~ invariant mass distributions after requiring \rriK± n o — 0.896| < 0.06 
GeV/c 2 are shown in Fig. H The spectra are fitted with Monte Carlo determined p 
shapes and Legendre background polynomials, as shown in Figs. a) and (b), and 
the fit yields 92 ± 19 events at y/s = 3.686 GeV and 5 ± 4 events at a/s = 3.65 GeV 
for K*(892) + Ii~p° ±c.c.. Similarly the 7r ± 7r° invariant mass distributions are shown 
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FIG. 4: K ± n T invariant mass distributions for K + candidate events from (a) 
s/s = 3.686 GeV data and (b) s/s = 3.65 GeV data. The squares with error bars are data, 
the histograms are the fit, and the dashed lines are the sum of the Kn mass distribution 
from 7T° sideband events and a Legendre polynomial. 



FIG. 5: K±ir° invariant mass distributions for K + K~7t + tt~tt 0 candidate events from (a) 
s/s = 3.686 GeV data and (b) s/s = 3.65 GeV data. The squares with error bars are data, 
the histograms are the fit, and the dashed lines are the sum of the Kn mass distribution 
from 7T° sideband events and a Legendre polynomial. 
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in Fig.ITJfor the y/s = 3.686 GeV and y/s = 3.65 GeV data samples. The fits, shown 
in Figs. Efa) and (b), yield 142 ± 23 events at y/s = 3.686 GeV and 6 ± 4 events 
at y/s = 3.65 GeV for K *(892)°A' _ p + + c.c.. After subtracting the continuum 
contributions, we obtain 78 ± 23 ip(2S) — > K*(892) + K~p° + c.c. events with a 
detection efficiency of (2.32 ± 0.05)% and 125 ± 26 -0(2 S) —> K* (892)° K~ p + + c.c 
events with a detection efficiency of (2.24 ± 0.05)%. 




u + u invariant mass (GeV/c 2 ) 


FIG. 6: 7r + 7r - invariant mass distributions from (a) /s = 3.686 GeV data and (b) 

y/s = 3.65 GeV data for J\*(892) + + c.c. candidate events. The squares with 
error bars are data, the histograms are the fit, and the dashed curves are background 
shapes from the fit. 

After requiring \m Kn — 0.896| < 0.06 GeV/c 2 and [m n7r — m p \ < 0.15 GeV/c 2 , the 
Dalitz plots of V>(2 S) -»• K*(892)+K~p° + c.c. and ijj(2S) -»• K*(892)°K~p+ + c.c. 
candidates are shown in Figs. H a) and (b), where no further clear structures are 
observed. 


V. SYSTEMATIC ERRORS 

Systematic errors on the branching fractions mainly originate from the MC statis¬ 
tics, the track error matrix, the kinematic fit, particle identification, the photon ef¬ 
ficiency, the uncertainty of the branching fractions of the intermediate states (from 
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ji + jt° or 7i 7i° invariant mass (GeV/c 2 ) 7i + n° or n n° invariant mass (GeV/c 2 ) 


FIG. 7: invariant mass distributions from (a) y/s = 3.686 GeV data and (b) 

y/s = 3.65 GeV data for A"*(892)°/\ _ 7r + 7r° + c.c. candidate events. The squares with 
error bars are data, the histograms are the fit, and the dashed curves are background 
shapes from the fit. 

PDG) Q], the K*{ 892) simulation, the fitting, and the total number of i^(2S) 
events. 

1. The MDC tracking efficiency was measured using clean channels like J/tj) —> 
AA and '0(25') —> 7r + 7r _ J/0, J/0 —> /r + /iV It is found that the MC simulation 
agrees with data within (1 — 2)% for each charged track. Therefore, 8% is taken 
as the systematic error for events with four charged tracks. 

2. The photon detection efficiency was studied using different methods with 
J/'tp —> 7T + 7T~7i° events js|, and the difference between data and MC simula¬ 
tion is about 2% for each photon. The systematic error due to the differences 
between data and Monte Carlo simulation of fake photons and the reconstruc¬ 
tion of the 7r° is less than 2.5% jnj. We take 5% as the systematic error for 
channels with two photons. 

3. The systematic error associated with the kinematic fit is caused by differences 
between the momenta and error matrices from track fitting of the charged 
tracks and the energies and the directions of the neutral tracks for data and 





























13 


> 

a) 

C3 


8 - 


6 - 


4 - 



CM 


• (°) 

0 8 

CD 

"O 

0 

CC 

•; ■ ■ ••• (b) 

*• * .. • * 

• n 


cr 

w 6 
0 

0 

0 

E 

* • • • t • 

* . • • • 

'o. 

o 

4 

* * • •• 


+ 

Cl 

o 

• • • v • * * * • •• 

_._._,_i_,_,_,_i_,_,_,_i_ 

■5K- 

2 

_1_._._1_._._._1_._._1_1_ 


K'p° or K + p° mass squared (GeV/c 2 ) 2 


4 6 8 

K"n nr K + n mass annarnri rGeV/n 2 ^ 2 


FIG. 8: Dalitz plots of (a) ip(2S) —► K*(892) + K p° + c.c. and (b) i/)(2S) 

K* (892)° K~ p + + c.c.. 


Monte Carlo data. For the channels analyzed, Prob(\ 2 i 4) > 0.01 is required, 
and we take 4% as the systematic error from the kinematic fit 12]. Fig. El 
shows the comparison of the y 2 distributions for data and MC sample, where 
we require |m 77 — 0.1351 < 0.03 GeV/c 2 and subtract tc° sideband background. 
MC simulation agrees with data within large statistical uncertainty of the data 
sample. 


4. The background uncertainties in the K + K~TT + 7r~n °, u>K + K~, cn/o(1710), 
Ji*(892)°/i _ 7T + 7r° + c.c., iC*(892) + iF - 7r + 7r - + c.c., K* (892) + K~p° + c.c., and 
K* (892)° K~p + + c.c. channels are estimated to be about 1.5%, 1.6%, 1.6%, 
7.1%, 9.5%, 11.4%, and 11.3%, respectively, by changing the order of the 
polynomial and the fitting range used. Varying the mass and width of the 
/o(1710) 0 in the fit yields a change in the number of cn/o(1710) events 
by 8.6%. Together with the background uncertainty, we take 8.8% as the 
systematic error for ^(2 S) —> u;/ 0 (1710) due to fitting. 


5. Pure 7r and K samples were selected, and the particle identification efficiency 
was measured as a function of the track momentum. On the average, a 
1.3% efficiency difference per n track and a 1.0% difference per K track are 
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FIG. 9: The \' 2 distribution from the kinematic fit. The squares with error bars are data, 
and the histogram is MC simulated ip(2S) —> K + A" _ 7r + 7r“7r°. 


observed between data and MC simulation. The total systematic error for 
K + K~ 7r + 7r _ 7r° is taken as 4.6%. 


6. The K*( 892) is simulated with a P-wave relativistic Breit-Wigner function, 

2 2 “ . ^ 

with the width T = T 0 — . , V? — , where r is the interaction radius and the 

u m 1 +r A p z [po J 7 

value (3.4 ± 0.6 ± 0.3) (GeV/c) _1 measured by a K~n + scattering experiment 


Q 


is used as an estimation of the interaction radius. Varying the value of r by 
la, the detection efficiencies for ij}(2S) decaying to A"*(892)°A' _ 7r + 7r° + c.c., 
iP*(892) + iP _ 7r + 7r _ + c.c., K* (892) + K~ p° + c.c., and K*(892)°K~p + + c.c. 
changed by 3.0% for the first two modes and 7.0% for the last two modes, 
which are taken as systematic errors for the uncertainty of the r value. 


7. For the JF*(892)°/v ~7r + 7r° + c.c. decay mode, there are backgrounds that 
also peak in the signal region of the Ki r mass plot. The largest is from 
7/l*(892)°A'~7r + + c.c. events combined with fake photons. We use n° side¬ 
band events to estimate these background contributions. In the analysis, the 
7T° sideband range is from 3 to 5 a, where a is 15 MeV/c 2 . If the n° sideband 
range is chosen between 5 and 7 cr, the difference with the standard sideband 
is 8.0%. If the 7r° sideband range is from 2 to 3 a, the difference is 13.7%. So 
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we take 13.7% as the systematic error due to the 7r° sideband definition in the 
fitting. Using a similar procedure for K*(892) + K + c.c. decay, 8.1% is 
determined as the systematic error. 

8. In calculating the Born order cross section for e + e“ — > K + K~ n + tt~ 7t° , the 
error for the integrated luminosity is 4%. The other systematic errors are 
similar to those for ip {2S) —> K + K~n + TT~n 0 . The total systematic error for 
a(e + e~ —> K + K~tt + tt~tt°) is 12.1%. 

Table □ lists all the systematic errors from different sources, the total system¬ 
atic errors for ip(2S) —► K + K~7t + 7t~7t 0 , ip{2S) —> uK + K~, ip(2S) —■> u;/ 0 (1710), 
ip(2S) -»• K*°K~tt + tt°+c.c., ip{2S) -»• K*+ K~tt + tt-+c.c., ip(2S) -»• K* + K~p°+c.c. 
and ip{2S) -»• K*°K~p + + c.c. are 12.1%, 12.2%, 15.1%, 19.9%, 17.7%, 18.1% and 
18.1% respectively. 


TABLE I: Summary of systematic errors (%). Errors common to all modes are only listed 
once. 


Source 

K + K 7T+7T 7T° 

wK+K- 

w/ 0 (1710) K*°K-ir+n° 

+C.C. 

K* + K-TT+Tt- 

+C.C. 

K*+K~p° 

+ C.C. 

K*°K~p+ 

+C.C. 

MC statistics 

1.4 

1.5 

2.3 

2.0 

1.8 

2.2 

2.3 

Tracking efficiency 




8.0 




Kinematic fit 




4.0 




PID efficiency 




4.6 




Photon ID and 








7r° reconstruction 




5.0 




Fitting 

1.5 

1.6 

00 

oo 

7.1 

9.5 

11.4 

11.3 

Branching fractions 


0.8 

0.8 





r uncertainty 




3.0 

3.0 

7.0 

7.0 

7T° sideband 




13.7 

8.1 



N ip(2S) 




4.0 




Sum 

12.1 

12.2 

15.1 

19.9 

17.7 

18.1 

18.1 


VI. RESULTS AND DISCUSSION 

To obtain the branching fraction of ip(2S) —> A", we must subtract the contri¬ 
bution from the continuum process. This is estimated using continuum data at 
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a/s = 3.65 GeV, normalized by /: 

r w _cont 

r _ *-3.686 x ^3.686 

J r w ^.cont ’ 

*-3.650 X CT 3 650 

where £^ is the integrated luminosity at a/s, and cx(/U is the Born order cross 
section of the continuum process at a/s, which is s dependent and can be expressed 
in terms of a form factor A(s): 

A'/r/y^ 

<r = ^ x i-ns)r, a) 

where a is the QED fine structure constant. Assuming |A(s)| oc 1/s, we get / = 
2.89. The branching fraction of il>(2S) —> A" can be calculated from 


B[i}{2S) -> X] 


jsrohs _ Arc 

iv 3.686 JV 3.650 


robs 


x/ 


^y(2S) 


X 


Artot 

iV b(2 S) 


X 


B( X 


Y)’ 


where X is the intermediate state and Y is the final state. 

Using the numbers obtained above and listed in Table El we get 


B{^{2S) -> K + K~n + n~n°) = (1.17 ± 0.10 ± 0.15) x 10~ 3 , 

B(${2S) -> uK + K~) = (2.38 ± 0.37 ± 0.29) x 10” 4 , 

B(i)(2S) -> w/ 0 (1710), / 0 (1710) -»■ K+K~) = (5.9 ± 2.0 ± 0.9) x 10" 5 , 

B{$(2S) -> A'*(892)°A'”7r + 7r° + c.c.) = (8.6 ± 1.3 ± 1.8) x 10“ 4 , 

B(il)(2S) -> A'*(892) + A'~7r + 7r _ + c.c.) = (9.6 ± 2.2 ± 1.7) x 10“ 4 , 

B(^{2S) K* (892) + K~ p° + c.c.) = (7.3 ± 2.2 ± 1.4) x 10“ 4 , 

B(^{2S) -> K* (892)° K~ p + + c.c.) = (6.1 ± 1.3 ± 1.2) x 10~ 4 , 

where the first errors are statistical and the second are systematic. The measured 
i)(2S) —> K + K~TT + 'K~n° branching fraction agrees well with the value of (1.27 ± 
0.05 ± 0.10) x 10 -3 obtained by CLEO j3], and the -0(25) —> uK + K~ branching 
fraction agrees with BESI Q and CLEO results within 1.5cr. The other five modes 
are first observations. 

The Born order cross section for e + e _ —> Ii + K~tt + 'k~ti g at a/s = 3.65 GeV is 

jyobs 

<7 e+e-^K+K-TT+-K-ir° = ^ = ± 35 ± 21) pb, 
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TABLE II: Numbers used in the branching fraction calculations. The number of events 
due to tp(2S) decay, N° b ^s) ■ is computed according to N§ b ^ S6 — f x No¬ 


Quantity 

K + K 7r"*”7T 7T° 

ujK+K- 

aJ /o(1710) 

K*°K~ TT+1T 0 

K*+K-tt+tt- 

K*+K~p° 

K* 0 K~p+ 





+C.C. 

+C.C. 

-\~c.c. 

+C.C. 

pjobs 

iV 3.65 

35 ± 7 

O+i.3 

U -0 

- 

15 ±5 

6 ± 5 

5 ± 4 

6 ± 4 

Mobs 

JV 3.686 

698 ± 41 

78 ± 11 

18.9 ±6.2 

281 ± 30 

150 ± 26 

92 ± 19 

142 ± 23 

lyobs 

>(2 S) 

597 ± 46 


18.9 ±6.2 

238 ± 34 

133 ± 30 

78 ± 23 

125 ± 26 

e{%) 

3.68 ±0.05 

2.66 ±0.04 

2.62 ±0.06 

3.00 ±0.06 

3.02 ±0.06 

2.32 ± 0.05 

2.24 ±0.05 

^(2S)(10 6 ) 




14(1 ± 4%) 





where e is the detection efficiency obtained from the MC simulation, C is the inte¬ 
grated luminosity and 1 + S is the radiative correction factor which is 1.23 0- 
To tost; the 12% rule, we ctlso list in Tctble nm the ratio Q h of the and 

J/ifj branching fractions for the three channels —> K + K~7r + Tr~n°, ip{2S) —> 

ujK+K~, and ip(2S) -*• u f 0 (1710), / 0 (1710) -»• K+K~. B{J/^ -»■ K + K~7r + 7r~7r°) 
is taken from the PDG jXQj 5 wliile the values on —» ujK + K~ ) and B(J/ip —> 


4(1710)1 come ,om Ref. 0, a recent mea_ t us ing a pa rti a> wave ana, ysis 
and the BESII J/ip sample. All three modes obey the 12% rule within la. The four 
modes with K* we observed in if)(2S) decays are not measured in J/-0 decays, so 
Qh values can not be computed. 


TABLE III: Branching fractions for ip(2S) and «//?/> hadronic decays and Qh values. The 
errors are the quadratic sum of the statistical and systematic errors. 


Channel 

4 ) 

4 ) 

Q h (%) 


11.7 ± 1.8 

120 ± 28 [10] 

9.8 ±2.8 

loK + K~ 

2.38 ± 0.47 

16.8 ±2.1 [13] 14.2 ±3.4 

w/ 0 (1710) -> uK+K~ 

0.59 ±0.22 

6.6 ± 1.3 [13] 

8.9 ± 3.8 
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